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Abstract 

ZnO  nanowires  with  the  [0110]  growth  orientation  undergo  a  reversible 
phase  transformation  from  wurtzite  (WZ)  to  a  graphitic  phase  (HX)  under 
tensile  loading,  leading  to  a  pseudoelastic  behavior  with  recoverable  strains 
up  to  16%.  Here,  we  report  that  this  phase  transform  causes  a  novel  transition 
in  thermal  response.  Molecular  dynamics  simulations  with  the  Green-Kubo 
approach  are  carried  out  to  determine  the  thermal  responses  of  wires  in  the 
18.95-40.81  A  size  range.  Results  obtained  show  that  the  thermal 
conductivity  in  the  unstressed  WZ  state  is  8. 3-8. 6  W  m-1  K-1  for  the  sizes 
considered  and  is  an  order  of  magnitude  lower  than  the  corresponding  bulk 
value.  Under  loading,  elastic  stretching  and  the  formation  of  interfaces  cause 
the  thermal  conductivity  to  first  decrease  significantly.  As  the  transformation 
progresses,  the  conductivity  increases  rapidly  to  10.1-11.1  W  m-1  K-1,  or 
20.5-28.5%  higher  than  that  of  the  initial  WZ-structured  wires.  The 
enhancement  is  primarily  due  to  an  increase  in  atomic  packing  density,  lower 
anharmonic  coupling  of  phonons  and  higher  surface  specularities  of  the 
HX-structured  wires.  This  phenomenon  offers  a  means  for  tuning  the  thermal 
behavior  of  ZnO  nanowires. 


1.  Introduction 

A  generation  of  semiconducting  nanostructures  such  as 
nanowires,  nanotubes,  nanobelts  and  nanorods  has  been 
developed  recently.  As  functional  building  blocks,  these 
nanocomponents  have  a  wide  variety  of  potential  applications 
as  piezoelectric  generators,  force  sensors,  resonators,  chemical 
and  biomolecular  sensors,  transparent  conductors,  catalysts 
and  nanoelectronic  and  nanophotonic  components  [1-6]. 
Challenges  exist  for  the  integration  of  these  structures  in 
nanoelectronic  systems  due  to  high  electric  field  strength  and 
high  heat  flux  at  the  nanoscale.  In  particular,  heat  dissipation 
and  thermomechanical  reliability  are  of  great  concern.  One 
issue  in  this  regard  is  that  the  thermal  conductivities  of 
nanostructures  are  1-2  orders  of  magnitude  lower  than 
corresponding  bulk  values  due  to  boundary  scattering  of 
phonons  [7-15]. 

The  coupling  between  the  mechanical  and  thermal 
responses  of  these  materials  provides  a  mechanism  for 
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tuning  and  increasing  the  thermal  conductivity  through  the 
application  of  mechanical  input.  Transitions  in  thermal 
conductivity  and  the  coefficient  of  thermal  expansion  under 
pressure  have  been  reported  for  the  quartz-coesite,  olivine- 
y -spinel,  coesite-stishovite,  B1-B2  and  pyroxene-garnet 
transformations  in  minerals  [16-19].  It  has  also  been 
shown  that  the  thermoelectric  properties  of  bulk  antimony 
bismuth  telluride  can  be  tuned  and  optimized  through  applied 
pressure  [20].  Although  thermal  transport  at  the  macro-  and 
nanoscales  is  well  characterized,  the  effect  of  stress  on  the 
thermal  response  of  nanostructures  has  not  been  extensively 
studied.  The  analysis  of  Picu  et  al  [21]  on  the  dependence  of 
thermal  conductivity  of  argon  on  hydrostatic  stress  provides 
some  insight,  without  accounting  for  the  effects  of  phase 
transformation  and  surfaces  which  play  a  dominant  role  at  the 
nanoscale  since  thermal  processes  are  dominated  by  surface 
scattering  of  phonons. 

Recently,  we  reported  a  novel  reversible  phase  transfor¬ 
mation  from  wurtzite  (WZ,  P6i,mc)  to  a  previously  unknown 
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graphitic  polymorph  (HX,  P^lmmc)  in  [01 10] -oriented  ZnO 
nanowires  under  tensile  loading  [22,  23].  This  discovery  has 
subsequently  been  confirmed  in  [0001] -oriented  nanoplates 
with  thicknesses  below  20  A  (cf  [24])  and  in  nanowires  with 
diameters  below  13  A  [25].  Here,  we  analyze  the  effect  of  this 
transformation  on  the  thermal  conductivity  of  such  nano  wires. 
The  analysis  focuses  on  variations  in  thermal  conductivity  dur¬ 
ing  the  elastic  stretching  of  the  initial  WZ- structured  wire,  the 
phase  transformation  into  the  HX  structure,  and  upon  comple¬ 
tion  of  the  phase  transformation.  The  characterization  of  the 
thermal  response  as  a  function  of  strain  here  is  important  since 
axial  elongation  is  one  of  the  most  relevant  modes  of  deforma¬ 
tion  for  slender  ID  nanomaterials. 

2.  Computational  framework 

Molecular  dynamics  (MD)  simulations  using  the  Buckingham 
potential  with  charge  interactions  [26,  27]  are  carried  out.  The 
Buckingham  potential  for  ZnO  has  been  shown  to  accurately 
predict  the  equilibrium  lattice  energy,  cell  parameters,  elastic 
and  dielectric  constants.  Extensive  perfect  lattice,  defect 
and  monovalent  ion  incorporation  simulations  have  been 
successfully  carried  out  using  this  potential  [26,  28,  29].  The 
potential  also  effectively  predicts  surface  properties  such  as 
surface  energies  [29].  This  aspect  is  especially  important  in 
the  simulations  for  nanowires  whose  high  surface-to-volume 
ratios  are  known  to  significantly  affect  behavior.  Recently,  MD 
simulations  on  nanowires  using  this  potential  yielded  thermal 
and  mechanical  behaviors  that  are  in  excellent  agreement  with 
the  predictions  of  first-principles  calculations  [14,  23]. 

The  nanowires  considered  are  single-crystalline  and 
wurtzite- structured,  with  lattice  constants  a  =  3.249  A  and 
c  =  5.206  A  and  a  growth  direction  along  the  [0110] 
axis  [4,  5].  Three  different  cross-sectional  sizes  (21.22  x  18.95, 
31.02  x  29.42  and  40.81  x  39.89  A2)  are  considered.  The 
smallest  cross-sectional  size  (21.22  A  x  18.95  A)  is  chosen 
such  that  the  short-range  cutoff  distance  in  the  Buckingham 
potential  [26,  27]  is  smaller  than  the  smallest  wire  dimension 
and  long-range  interactions  are  properly  considered  [30]. 
Periodic  boundary  conditions  are  specified  in  the  axial 
direction.  Calculations  with  different  computational  cell  sizes 
show  that  any  length  greater  than  100  A,  irrespective  of  the 
cross-section  size,  is  sufficient  to  avoid  image  effects  [14,  15]. 
Here,  a  periodic  computational  cell  length  of  150.83  A  is  used 
for  all  the  cross  sections  analyzed. 

The  deformation  of  the  nanowires  is  approximated  as  a 
quasi-static  process  with  each  deformation  increment  achieved 
though  successive  loading  and  equilibration  steps  using  a 
combination  of  algorithms  for  NPT  and  NVE  ensembles  [31]. 
Specifically  in  each  deformation  increment,  stretching  at  a 
specified  rate  of  0.005  ps-1  is  first  carried  out  for  0.5  ps  using 
a  modified  version  of  the  NPT  algorithm  of  Melchionna  et  al 
[32,  33].  Subsequently,  with  the  strain  maintained  constant, 
the  nanowire  is  relaxed  for  3  ps  via  an  algorithm  for  NVE 
ensemble  [31]  at  the  specified  temperature.  This  equilibration 
duration  is  chosen  such  that  a  statistically  steady  state  is 
reached  and  no  further  structural  changes  occur. 

Since  the  loading  proceeds  in  a  series  of  equilibration 
steps,  an  equilibrium  approach  such  as  the  Green-Kubo 
method  can  be  used  to  determine  the  thermal  response  of 


Figure  1.  Stress  and  thermal  conductivity  as  functions  of  applied 
strain  for  a  31.02  A  x  29.42  A  under  tensile  loading  at  500  K. 

(This  figure  is  in  colour  only  in  the  electronic  version) 

the  nanowires  as  a  function  of  strain.  Also,  since  the  MD 
calculations  for  determining  heat  flux  do  not  require  any  a 
priori  knowledge  of  the  heat  transport  processes  or  material 
configuration,  this  approach  is  ideal  for  evaluating  heat  transfer 
in  materials  undergoing  deformation  or  phase  transformation. 
At  each  strain,  the  thermal  conductivity  can  be  calculated 
according  to  the  fluctuation-dissipation  theorem  as 

Kw  =  777 ^  f  WW0)>  dt,  (1) 

F/C  Jo 

where  V  is  system  volume,  T  is  temperature,  k b  is  the 
Boltzmann  constant,  J^{t)  is  the  /xth  (/x,  v  =  1,2,3) 

component  of  the  heat  current  and  0))  is  the 

autocorrelation  function  for  J^it)  with  ( )  denoting  ensemble 
time  average.  The  analysis  is  conducted  at  500  K  which 
is  above  the  Debye  temperature  0D  (=420  K  for  ZnO). 
Consequently,  temperature  can  be  calculated  from 

N 

\NkBT  =  (2) 

i=  1 

with  N  being  the  number  of  atoms  in  the  system  and  w,- 
and  Vi  being,  respectively,  the  mass  and  velocity  of  atom 
i.  Temperatures  below  $D,  for  which  quantum  mechanical 
corrections  may  be  needed,  are  not  considered  here. 

3.  Results  and  discussion 

3.1.  Variation  of  thermal  conductivity  with  applied  loading 

Figure  1  shows  stress  and  thermal  conductivity  as  functions  of 
strain  for  a  31.02  A  x  29.42  A  wire.  Three  distinct  stages  of 
deformation  can  be  identified  from  the  stress-strain  response: 
(1)  elastic  stretching  of  the  defect-free  WZ- structured  wire  (A- 
B),  (2)  transformation  from  WZ  to  HX  (B-C)  and  (3)  elastic 
stretching  of  the  defect-free  HX- structured  wire  (C-D)  after 
completion  of  the  transformation  [22,  23]. 

3.1.1.  WZ-structured  wire  prior  to  initiation  of  phase 
transformation.  The  thermal  conductivity  (/c)  of  a  unstressed 
31.02  A  x  29.42  A  nanowire  at  500  K  is  8.52  Wm"1  K_1, 
an  order  of  magnitude  lower  than  that  for  bulk  ZnO 
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WZ  Interface  HX 


Figure  3.  Arrangement  of  atoms  on  a  (2110)  plane  across  an 
interfacial  region. 


Figure  2.  WZ,  HX  and  interface  regions  in  a  31.02  A  x  29.42  A 
wire  at  (a)  e  =  6.5%,  (b)  £  =  8.3%  and  (c)  £  =  10.5%. 


(~100  W  m_1  K_1)  [34].  This  significant  difference  is 
associated  with  significant  boundary  scattering  of  phonons 
due  to  the  high  surface-to-volume  ratios  of  the  nanowires 
and  changes  in  phonon  spectrum  for  the  quasi- ID  wire 
structures.  Specifically,  the  relatively  large  fractions  of  surface 
atoms  enhance  surface  scattering  of  phonons  and  decrease 
the  phonon  mean  free  path,  resulting  in  lower  conductivity 
which  is  proportional  to  the  mean  free  path.  In  contrast,  the 
characteristic  length  for  bulk  materials  is  much  longer  and  the 
effect  of  boundary  scattering  is  negligible,  resulting  in  longer 
mean  free  paths  and  much  higher  conductivity  values  [14,  15]. 
Furthermore,  changes  in  the  phonon  spectrum  at  the  nanoscale 
modify  phonon  group  velocities  and  scattering  mechanisms 
and  can  also  lead  to  lower  thermal  conductivity  values  and  size 
effects  [7,  35-37]. 

As  deformation  progresses,  k  decreases  from  its  initial 
value  of  8.52  W  m_1  K-1  by  19.6%  as  the  strain  increases  to 
0.049.  During  this  stage,  the  nanowire  is  fully  within  the  WZ 
structure  and  the  decrease  in  k  is  primarily  due  to  the  nonlinear 
elastic  behavior  of  WZ  along  the  [0110]  orientation,  although 
the  increase  in  lattice  anharmonicity  at  finite  values  of  strain 
also  has  a  contribution  [21].  The  nonlinear  elastic  response 
can  be  quantified  through  the  rate  of  change  of  the  uniaxial 
modulus 


d2a  _  1  d 3U 
d s2  ~  V  de5"’ 


(3) 


where  V  is  the  initial  volume  of  the  wire,  U  is  the  strain  energy, 
and  a  and  e  are,  respectively,  the  stress  and  strain  along  the 
wire.  For  the  [0110]  nanowire  analyzed,  E '  =  —1428.7  GPa 
and  the  modulus  E  decreases  by  35%  up  to  the  strain  at 
transformation  initiation  (e  =  0.049).  The  elastic  stretching 
also  causes  the  volume  of  the  wire  to  increase  and  the  mass 
density  p  to  decrease.  Both  changes  combine  to  cause  the 
stress  wave  speed  or  the  phonon  group  velocity  v  (=y/E/p)  to 
decrease  by  18.9%.  Additionally,  the  decrease  in  density  along 
with  the  lattice  distortion  increases  lattice  anharmonicity  and, 
consequently,  causes  the  mean  free  path  (A)  of  the  phonons  to 
decrease.  At  temperatures  above  the  Debye  temperature,  the 
kinetic  theory  of  fluids  relates  thermal  conductivity  to  phonon 
mean  free  path  through  [38] 


k  =  \Cvv  A, 


(4) 


where  Cv  is  specific  heat.  Obviously,  the  conductivity 
decreases  as  the  WZ-structured  nanowire  is  stretched. 


Figure  4.  Fractions  of  atoms  in  WZ,  HX  and  interface  regions  as 
functions  of  strain  for  a  31.02  A  x  29.42  A  wire. 


3.1.2.  WZ/HX- structured  wire  during  phase  transformation. 
Upon  the  initiation  of  phase  transformation  at  s  =  0.049,  the 
stress  drops  precipitously,  reflecting  a  relaxation  of  the  wire 
structure.  Figure  2  shows  the  configurations  of  the  wire  at 
s  =  0.065,  0.083  and  0.105.  Three  zones,  corresponding  to 
WZ,  HX  and  a  transitional  interface  between  the  two  phases, 
are  seen.  The  nucleation  of  HX  is  at  the  surfaces.  The 
atomic  structure  on  the  (2110)  plane  in  the  interfacial  region 
is  shown  in  figure  3.  Obviously,  the  interface  is  coherent 
and  the  orientation  relationships  are  maintained  in  both  phases 
across  the  interface.  The  fractions  of  atoms  in  the  HX,  WZ  and 
interfacial  regions  are  shown  in  figure  4  as  functions  of  strain. 
Initially  (0.049  <  s  <  0.072,  strain  range  from  transformation 
initiation  to  the  formation  of  the  largest  interface  fraction),  the 
fraction  of  the  HX  phase  is  lower  than  that  of  the  WZ  phase  and 
the  thermal  response  is  dominated  by  the  conductivity  of  the 
WZ  phase  and  the  thermal  resistance  of  the  interface.  In  this 
stage,  k  decreases  by  18.7%  from  6.86  to  5.58  W  m-1  K-1  as 
the  strain  increases  from  0.049  to  0.072.  This  decrease  is  due 
to  the  low  thermal  transmission  coefficient  of  the  interface  and 
enhanced  phonon  scattering  due  to  (1)  discontinuities  in  the 
lattice  across  the  interface  and  (2)  lattice  distortion  associated 
with  the  elastic  strain  field  near  the  interface  [39,  40].  Note 
that  over  this  strain  interval  the  fraction  of  interfacial  atoms 
increases  from  0. 12  to  a  maximum  of  0. 19.  Since  the  interfacial 
thermal  resistance  is  proportional  to  the  area  of  interfaces 
(quantified  here  through  the  fraction  of  interfacial  atoms),  k 
decreases  as  strain  increases.  Continuation  of  loading  beyond 
s  =  0.072  causes  the  fraction  of  HX  to  exceed  that  of 
WZ  (figures  2(b),  (c)  and  4),  resulting  in  the  overall  thermal 
response  of  the  wire  to  be  dominated  by  the  HX  phase.  At 
the  completion  of  the  WZ  -a-  HX  transformation  (s  =  0.113), 
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the  conductivity  of  the  HX- structured  wire  is  9.4  W  m_1  K-1 
or  68.3%  higher  than  the  lowest  value  of  5.58  W  m-1  K-1 
at  e  =  0.072.  This  increase  in  k  is  due  to  several  factors, 
including  (1)  the  higher  thermal  conductivity  of  the  HX  phase 
(details  later),  (2)  significant  decrease  (43.3%)  in  the  fraction 
of  interface  atoms  as  the  strain  increases  from  0.072  to  0.113 
and  (3)  the  alignment  of  the  interfaces  along  the  direction  of 
heat  flow  (wire  axis)  as  deformation  progresses  (figures  2(b) 
and  (c)).  Arising  primarily  from  bonding  state  non-uniformity 
and  lattice  strains,  interfacial  thermal  resistance  is  expected  to 
be  much  lower  for  directions  parallel  to  an  interface  than  that 
for  the  direction  perpendicular  to  it. 


Figure  5.  Arrangement  of  atoms  in  the  interior  and  on  the  surfaces  of 
a  31.02  A  x  29.42  A  wire,  the  images  correspond  to  (a)  bulk  WZ, 

(b)  WZ-structured  wire  with  LY-structured  surfaces  after  initial 
relaxation  and  (c)  HX- structured  wire  after  completion  of  the 
WZ  — >  HX  transformation. 


3.1.3.  HX-structured  wire  after  transformation.  The  WZ  —> 
HX  transformation  completes  at  a  strain  of  e  =  0.113.  At 
this  stage,  some  point  defects  exist  due  to  local  variations  in 
the  lattices.  Further  loading  causes  the  defects  to  disappear, 
leading  to  a  defect-free  HX  structure  at  s  =  0.125.  The  value 
of  thermal  conductivity  obtained  for  this  defect-free  structure  is 
reported  as  the  characteristic  thermal  property  of  the  HX  wires. 
Further  loading  causes  the  conductivity  to  decrease,  reflecting 
the  nonlinear  elastic  behavior  of  the  HX-structured  wires  (not 
shown).  The  thermal  conductivity  of  the  HX-structured  wire 
in  this  state  is  10.7  W  m-1  K-1,  which  is  25.6%  higher  than 
that  of  the  unstressed  WZ  wire.  Such  a  significant  increase 
in  k  can  be  qualitatively  explained  using  a  simplified  model 
for  thermal  conductivity  [18,  41].  At  temperatures  above  the 
Debye  temperature,  the  lattice  thermal  conductivity  is  limited 
by  anharmonic  coupling  of  phonons  and  can  be  expressed  as 
(cf  [18,  41]) 

rv3 

<5) 

where  r  is  the  average  interatomic  spacing  such  that  the  unit 
cell  volume  Vceii  =  r3,  v  is  the  average  wave  speed,  p  is  mass 
density,  y  is  the  thermal  Gruneisen  parameter,  T  is  temperature 
and  A  is  a  constant  representing  the  contributions  of  changes 
in  crystal  structure  (coordination,  bond  length,  etc.)  to  thermal 
conductivity.  During  the  WZ  -a-  HX  transformation,  the 
number  of  atoms  per  unit  cell  remains  the  same  and  the  average 
bond  length  increases  only  by  ~6%.  Hence,  the  change  in 
A  is  negligible  and,  at  a  given  temperature,  the  effect  of  the 
phase  transformation  on  thermal  conductivity  can  be  expressed 
through  the  change  in  density,  anharmonicity  of  the  lattice  and 
the  wave  speed  as 


The  wave  speed  v,  calculated  from  the  elastic  moduli 
and  densities  of  each  of  the  two  phases,  is  lower  for 
HX  than  for  WZ.  Therefore,  the  change  in  wave  speed 
tends  to  lower  thermal  conductivity  and  this  effect  is 
only  secondary.  Dominant  effects  that  cause  the  thermal 
conductivity  to  increase  come  from  the  higher  density  and 
higher  Gruneisen  parameter  of  HX.  During  the  transformation, 
the  lattice  parameter  along  the  [0001]  direction  (c)  decreases 
considerably,  causing  the  cell  volume  to  decrease  by  11.98% 
and  the  density  of  the  structure  to  increase  by  13.5%.  The 
efficient  packing  of  the  HX  structure  and  its  higher  density  lead 
to  a  lower  anharmonicity  (lower  thermal  Gruneisen  parameter) 


and,  ultimately,  the  higher  thermal  conductivity  [17].  For  the 
31.02  A  x  29.42  A  wire,  the  25.6%  increase  in  k  implies 
a  27.5%  decrease  in  y  relative  to  the  value  of  0.69  for 
WZ  at  500  K  [42].  This  decrease  is  consistent  with  the 
increase  in  packing  efficiency  of  the  HX  structure  over  the 
parent  WZ  structure.  A  quantitative  prediction  of  y  is  not 
feasible  here  since  the  model  is  strictly  valid  for  bulk  materials 
where  the  effect  of  surface  scattering  of  phonons  is  negligible. 
Nevertheless,  this  model  provides  a  qualitative  explanation  for 
the  increase  in  conductivity  observed  here. 

The  difference  in  surface  configurations  between  WZ- 
and  HX-structured  nanowires  also  contributes  to  the  observed 
increase  in  thermal  conductivity.  To  illustrate  this  issue, 
figures  5(a)  and  (b)  show  the  positions  of  atoms  on  layers 
perpendicular  to  the  [0001]  direction  for  bulk  WZ  and  a  WZ- 
structured  wire  and  figure  5(c)  shows  the  configuration  of  a 
HX-structured  nanowire.  For  WZ  wires,  the  imbalance  of  ionic 
forces  on  the  surfaces  due  to  the  reduced  number  of  neighbors 
and  surface  polarity  cause  extensive  surface  reconstruction 
relative  to  bulk  WZ.  Obviously,  this  reconstruction  entails  the 
contraction  of  surface  layers  and  merging  of  Zn  and  O  basal 
planes  (figures  5(a)  and  (b)),  resulting  in  a  layered  surface 
structure  (LY)  which  is  crystallographically  similar  to  the 
HX  structure  [43,  44].  Such  reconstructions  modify  atomic 
arrangement  on  surfaces  relative  to  the  core  of  the  nanowire 
and  alter  the  surface  scattering  behavior  of  phonons,  causing 
the  surface  specularity  and,  hence,  the  thermal  conductivity 
for  WZ-structured  nano  wire  to  be  much  lower  than  that  for 
bulk  WZ.  In  the  HX-structured  wire  in  figure  5(c),  the  atomic 
arrangement  on  surfaces  is  similar  to  that  in  the  core  and  hence 
the  surface  disorder  is  significantly  lower  than  that  of  the  WZ 
wire  in  figure  5(b).  As  a  result,  the  surface  specularity  for  HX 
wires  is  much  higher  than  that  of  WZ-structured  nanowires, 
resulting  in  much  lower  effects  of  boundary  scattering  and 
higher  conductivity  values.  The  effect  of  surface  specularity 
on  conductivity  can  be  characterized  through  equation  (4). 
Under  conditions  for  which  phonon-phonon  and  phonon- 
defect  interactions  are  negligible,  the  phonon  mean  free  path 
can  be  expressed  as  [45] 


where  p  is  the  probability  of  specular  scattering  (a  function 
of  surface  roughness  and  temperature)  and  d  is  the  effective 
size  of  the  nano  wire.  For  a  body  with  perfectly  specular 
(atomistically  smooth)  surfaces,  p  =  1.  Consequently,  the 
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Figure  6.  Thermal  conductivity  of  unstressed  WZ  and  transformed 
HX  wires  of  three  different  sizes. 


phonon  mean  free  path  is  infinite  and  the  thermal  conductivity 
is  not  affected  by  boundary  scattering.  On  the  other  hand, 
for  a  body  with  perfectly  diffuse  surfaces,  p  =  0  and  the 
mean  free  path  expression  reduces  to  the  Casimir  limit  such 
that  the  mean  free  path  of  the  phonons  is  limited  by  the  size 
of  the  nano  wire.  In  reality,  the  specularity  value  is  between 
the  limiting  cases  (0  <  p  <  1)  and  depends  significantly 
on  surface  attributes.  As  discussed  previously,  the  surface 
specularity  of  WZ  wires  is  much  lower  than  that  of  HX  wires. 
Consequently,  for  nanowires  with  the  same  characteristic  sizes, 
the  phonon  mean  free  path  and  therefore  the  conductivity  of  the 
HX  wires  is  higher  than  that  of  the  WZ  wires. 

3.2.  Size  dependence  of  thermal  conductivity 

The  bulk  and  surface  effects  discussed  above  combine  to 
give  rise  to  significant  dependence  of  behavior  on  wire  size. 
Figure  6  shows  the  conductivity  values  of  unstressed  WZ  and 
transformed  HX  wires  of  three  different  sizes.  For  WZ  wires, 
the  thermal  conductivity  is  in  the  range  of  8.3-8. 6  W  m-1  K-1 
which  is  an  order  of  magnitude  lower  than  that  for  bulk  ZnO 
and  decreases  by  3%  as  the  lateral  size  decreases  from  40.81  x 
39.89  to  21.22  x  18.95  A2.  This  trend  results  directly  from 
the  higher  surface-to-volume  ratios  and  smaller  mean  free  path 
at  the  smaller  sizes.  A  clear  dependence  of  conductivity  on 
size  is  also  seen  for  the  HX- structured  wires.  The  WZ-to-HX 
transformation  causes  the  conductivity  to  increase  to  10.1,  10.7 
and  11.1  Wm-1  K_1  for  the  21.22  x  18.95,  31.02  x  29.42 
and  40.81  x  39.89  A2  nanowires,  respectively.  The  20-28% 
increase  over  the  values  for  WZ  wires  is  associated  with  the 
lower  anharmonicity  due  to  the  higher  surface  specularity  of 
the  HX  wires. 

4.  Conclusions 

MD  simulations  with  the  Green-Kubo  approach  are  carried 
out  to  determine  the  thermomechanical  response  of  [0110]- 
oriented  ZnO  nanowires  under  tensile  loading.  Three  distinct 
stages  in  the  thermal  and  mechanical  responses  are  observed. 
The  mechanical  response  consists  of  (1)  elastic  stretching  of 
WZ,  (2)  phase  transformation  from  WZ  to  HX  and  (3)  elastic 
stretching  of  HX.  The  thermal  conductivity  of  the  nanowires 
is  a  function  of  applied  strain  in  each  of  the  three  stages. 


Over  the  size  range  of  18.95-40.81  A,  the  elastic  stretching 
in  the  WZ  structure  is  accompanied  by  19-27%  decrease  in  the 
thermal  conductivity.  The  nonlinearity  of  the  elastic  response 
of  the  wires  with  lower  stiffness  (hence  lower  wave  speeds) 
at  higher  strains  is  the  origin  of  this  effect.  The  formation  of 
interfaces  during  the  WZ-to-HX  phase  transformation  causes 
the  thermal  conductivity  to  decrease  by  another  15-19% 
as  the  transformation  progresses.  Upon  completion  of  the 
transformation,  the  thermal  conductivity  for  the  HX- structured 
wires  are  20.5-28.5%  higher  than  those  of  the  corresponding 
WZ- structured  wires.  This  increase  in  conductivity  is  due  to 
the  higher  atomic  packing  density,  lower  anharmonic  coupling 
of  phonons  and  higher  surface  specularity  of  the  HX  wires. 

Novel  phase  transformations  at  the  nanoscale,  like 
the  one  analyzed  here,  alter  thermal  conductivity  and 
provide  a  mechanism  for  controlling  or  ‘tuning’  the  thermal 
response  of  nanocomponents  through  the  application  of 
a  mechanical  input.  The  present  analysis  considers 
the  thermomechanical  response  of  ZnO  since  it  is  a 
versatile  semiconducting  material  with  applications  in  power 
generation,  sensing,  environmental  monitoring,  biomedical 
systems  and  communications  technology.  It  should  be  noted 
that  such  phase  transformations  can  also  occur  in  other  group 
IV,  III-V  and  II- VI  materials  including  InN,  GaN  and  SiC  [46], 
leading  to  a  similar  potentially  exciting  tunability  of  thermal 
responses  in  these  materials  as  well.  The  thermomechanical 
coupling  in  these  materials  may  therefore  provide  both 
opportunities  for  developing  ‘tunable’  functional  nanodevices 
and  challenges  for  ensuring  thermomechanical  reliability  and 
functionality  of  a  range  of  NEMS. 
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